NEW EQUATIONS FOR VACUUM CORRELATORS 

IN GAUGE THEORIES 



D.V.Antonov 

! Institute of Theoretical and Experimental Physics, 

Gn ■ B.Cheremushkinskaya, 25, 117218, Moscow, Russia 

t-H ■ 



> 

in 



Abstract 

Stochastic quantization [1] is applied to derivation of equations, connecting 
multilocal gauge-invariant correlators in different field theories. They include 
Abelian Higgs Model, QCD with spinless quarks at T = and T > and 
O ; QED, where spin effects are taken into account exactly. 
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One of the main problems of the modern quantum field theory is de- 
termination of nonperturbative effects in the models with nontrivial vacuum 
structure, like QCD. To investigate this problem the Method of Vacuum Cor- 
relators (MVC) [2-8, 17] (for a review see [5]) was suggested. In the framework 
of this method one, using Feynman-Schwinger path integral representation, 
extracts from quark- antiquark Green function all the dependence from 
vacuum gluonic fields in the form of an averaged Wilson loop. After that, 
through nonabelian Stokes theorem [3] and cumulant expansion [3,9] the ave- 
raged Wilson loop is expressed via an infinite set of irreducible gauge- 
invariant vacuum correlators (cumulants). Lattice data suggest, that the 
ensemble of gluonic fluctuations is predominantly Gaussian [5] . That is why, 
in order to calculate observable quantities (like confining potential [4,5], mass 
spectra of mesons, baryons and glueballs [5-7], Regge trajectories [7,8], tem- 
perature of deconfinement, meson and glueball screening masses at finite 
temperature [17], etc.) one may with a good accuracy neglect all the cumu- 
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lants higher than quadratic and to parametrize the latter by two independent 
Kronecker structures, where the space-time dependence of the corresponding 
coefficient functions is determined from the lattice data. 

However, the nonperturbative stochastic Euclidean QCD vacuum is de- 
scribed by the full set of cumulants, and, therefore, for the completeness 
of the theory, one needs the equations (derived from the Lagrangian), from 
which the cumulants may be obtained. Two alternative methods of deriva- 
tion of such equations, based on stochastic quantization [1], were suggested 
in [10]. First of them, exploiting Feynman-Schwinger path integral represen- 
tation, where Langevin time is considered as a proper time, was applied to 
the (p^ theory. However, it is inconvenient for gluodynamics, since it breaks 
down gauge invariance. That is why, in the latter case an alternative ap- 
proach, leading to integral-differential Volterra type-II equations, was used. 
In both cases it was shown, that the perturbative expansion of the obtained 
equations in the lowest orders reproduced correctly known results of stan- 
dard Feynman diagrammatic technique. An important point is that in the 
physical limit, when Langevin time t tends to infinity, the solutions of these 
equations yield exact correlators of Heisenberg field operators, averaged over 
the physical vacuum, whose structure may be complicated and unknown, e.g. 
in QCD. 

In this paper we demonstrate, how these two approaches work simul- 
taneously in gauge theories with matter fields, presenting the methods of 
derivation of equations for gauge-invariant correlators in different models. 
The sketch of the paper is the following. In section 2 we start with the 
Abelian Higgs Model (AHM), where classical solutions (Abrikosov-Nielsen- 
Olesen (ANO) strings) are present. These objects may be naturally included 
into obtained equations as initial conditions in the Cauchy problem for the 
Langevin equation, and, thus, we obtain a new method of quantization of 
classical solutions. This method allows one to derive correlators, containing 
arbitrary number of quantum fluctuations around classical solutions in any 
theory. We hope, that it will be especially useful in the case of gauge theories, 
since the obtained equations are explicitly gauge- invariant. Computation of 
quantum corrections to an instanton in the framework of this method will 
be a topic of a separate publication. Another point is that the suggested 
approach is the first attempt to gauge-invariant quantization of theories with 
spontaneous symmetry breaking and may be applied to quantization of more 
complicated theories, like Standard Model. In section 3 we exploit methods 
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of stochastic quantization of fermionic theories to derive equations for gauge- 
invariant correlators in QCD with spinless quarks. Then we demonstrate, 
how the obtained equations may be generalized to the finite temperature 
case. Finally, in section 4 we use the path integral approach to calculation 
of the propagator of a spinning particle in the Abelian background gauge 
field [11,12] to obtain equations in QED, where all spin effects are taken into 
account exactly. In Conclusion we summarize the main results of the paper 
and outline possible future developments. In the Appendix we obtain a Green 
function of the Langevin equation for F fIU (x : t) in the AHM. 



2 Abelian Higgs Model 

In this section we consider AHM, where exists a vacuum of classical solu 
tions (ANO strings) [13]. The Euclidean action of the AHM has the form 
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S = J dx 



(1) 



where D^ip = (d M — ieA^)p . Let us look for a solution of equations of motion 
in the cylindrical coordinates in the form A(r) = Aq{t), A4 = 0, <p =| <p \ e m0 , 
where | (p | is a function of r only. Then it is easy to check, that such a solution 
(ANO string along z-axis), corresponding to an n-quantum Abrikosov vortex 
of Ginzburg-Landau model, has the form 

Mr >-\£t(-i + ip-plnp), t«r«6, (2) 



where 5 = — is a Londons' depth of penetration of the field A u , f = — 

me " " A" => m 

is a correlation radius of fluctuations of the Higgs field, p = j,K\ is the 
Macdonald function. 

The Langevin equations, following from the action (1), have the form 

p = D 2 p — m 2 p + 2\p*p 2 + 77, (3) 

tp* = (D 2 ip)* - my + 2\iptp* 2 + 77*, (4) 

A li = d v F vli +j ti + r} fl , (5) 
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where = —ie((p*d^(p — (d^(p*)(p — 2ieA fl tp*tp) is a gauge-invariant current, 

< r](x,t)ri*(x',t') >= 28(x - x')8{t - t'), 

< rj^x, t)r] v {x\ t') >= 25^5(x - x')8{t - t'), (6) 

and < ... > means the averaging over both Gaussian stochastic ensembles of 
(77, 77*) and Tj^. 

Passing from equation (5) to the Langevin equation for F^x^t) and 
solving it via Fourier transformation with the initial condition A^(x,§) = 
= A^x), where A r {x) = A z (x) = A±{x) = and Ag(x) is defined through 
the formula (2), one obtains: 

t 

F liU (x,t) = j dt' J dyU^(y-x,t-t')(jp(y,t')+r)p(y,t'))+ 



+ / dy U^(y - x,t)A (y). 



(7) 



where 



48tt 2 



(Sua 8v0 — &vot 8np) 



+ 



2 / 1 2 \ 1 

+^2(^2/^ - ^y v )i^e~^ - 16n 2 5(y)j + ^[5^y u y a + 



i^yp 



(8) 



The details of derivation of the Green function (8) are presented in the Ap- 
pendix . 

Using Langevin time as a proper time and solving equation (3) via 
Feynman-Schwinger path integral representation with the initial condition 

<p(x,Q) = <po(x), I (fo \ 2 = we have: 



1 

<p(x,t) = J dt' J dy(Dz) xy K z (t,t')F z (t,t!)<S>(x,y)ri{y,t'), 



(9) 



where 



n d 4 z(n) 

(Dz) xy = lim n 77^, ,Ne = t-t', z^ = t')=y, z(£ = t)=x, 




t 



T z {t,t') = exp 2\J d£\ <p(z(£),£) 



v 



x 



$(x,y) = exp ie J d z ^ A ^(z 

y 



and we omitted in (9) the term / dy(Dz) xy K z (t,0)J r z (t,0)&(x,y)(po(y) , be- 
cause we shall see, that in the correlators (11), (12) (and all the other gauge- 
invariant quantities) this term produces non-gauge-invariant part of the cor- 
relator, which in the physical limit, t — > + oo, tends to zero. 
In the same way one may solve equation (4). 

To obtain an infinite system of exact equations for gauge-invariant 
correlators, containing the fields <p and </?*, one should multiply 
the right hand sides of equation (9) and the equation for <p* by 
r](x,t), r]*(x,t), r^(x,f), ip(x,i), <p*(x,i), j^(x,i), F^(x, t) as many times as 
needed, introducing corresponding parallel transporters in order to preserve 
gauge invariance, and average both sides of the obtained equations. The right 
hand sides of these equations may be most easily expressed through their left 
hand sides by introducing the generating functional 



where A: is a total number of fields cp and p* in the correlator, Wi{Cj) is an 
Abelian Wilson loop, and the form of a contour d is clarified below. After 
that one should apply to this generating functional cumulant expansion [3,9]. 

For example, let us consider the simplest case k = 2. This is the so-called 
bilocal approximation [10], when one assumes, that the ensemble of the fields 
p, (p*,jfj,, F^ v and stochastic noise fields is Gaussian, so that all the cumulants, 
higher than quadratic, are equal to zero. In this case we obtain 



Z= ]J < Fzi(ti,lfi)Wi(Ci)exp / dudt(I li (u,t)rj li (u,t)+ 



i=l 




(10) 



< p>*(x,i)&(x,x,r)ip(x,t) > 



5 



min(t,t) 

= 2 / dt'J dy{Dz) xy {Dz) xy K z {t,t , )K z {t,t') < T z (t,t')T- z (f,t')W yxxy >, 
o 

(11) 

< rf(x,l!)$(x,x,T)ip(x,t) >= - < ip*(x,t)$(x,x,T)r)(x,1?) >= 

= 2 / (Dz) xx K z (t, t') < T z {t, t')W xxx >, (12) 

where <$>(x,x,t) is a parallel transporter between the points x and x, given 
at the moment r of Langevin time (in particular, one may consider the case, 
when x = x and $>(x, x, r) = 1). Here 



Wyxxy - exp < 



le 



{|- x x 
ie J du^A^u.r) + j dz^ A ^{z 



(13) 
(14) 



Formulae (13) and (14) clarify the construction of the contours C\ and C 2 
in bilocal generating functional. Applying to it cumulant expansion, varying 
over M or N in the case of equations (12), putting 1^ = = = M = 
N = and rewriting W yxxy and W xxx through Stokes theorem, one obtains 
in the bilocal approximation: 

r 1 

< T z (t,t')T- z (t,t')W yxxy >= exp 2A / d£ <| | 2 > + 



+2A/de<l^(0,f) i 2 > + 



Si 



Si 



< T z {tJ)W xxx >= exp 



2\jdi<\ V (z(0,0 | 2 > + 



+ \ / ^iiW I d<r nn (u 2 ) < F Hl , I (u l ,t l )F mU2 (u2,t 2 ) > 



(15) 



(16) 



where 5i and S 2 are arbitrary surfaces, bounded by the contours C\ and C 2 
respectively, t\ and t 2 are some moments of Langevin time. 
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The solution (7) of the Langevin equation for F p v(x,t) yields three more 
equations for gauge-invariant correlators: 

t t 

< F fiV (x,i)F Xp (x,t) >= J dt' J d¥ J dydy U^ >a (y-x, t-t')U Xp ^{y -x, t-?)- 



o o 



< ja(y,t')j(3{yJ) > + < ja(y,t')r]p(yj) > + < jp(yj)va{y,t') > + 



+2S afi (y - yW ~ i') 



+ 



+ J dydyU^ a (y - x, t)U Xp ,p(y - x, t)A a (y)Ap(y), (17) 



i 

< Fp„(x, t)j\(x, t)>= J dt' J dy U^ a (y -x,t- t') 



+ < j\{%,t)r] a {y,t') > 



< F pi/ (x, t)rix(x, t) >= J dt' j dy U^ a (y -x,t- t') 



+ 28 aX 8(y - x)8(t' - t) 



< ja(y,t')j X (x,i) > + 

(18) 

< 3a(y,t')r]x(x,t) > + 

(19) 



In order to obtain a complete system of equations, let us pass from the 
Hamilton gauge to the Schwinger gauge, A p (x,t)(x — xq) p = (where xq is 
an arbitrary point), in which it is easy to express explicitly A p (x,t) through 
F pv {x,t) and to obtain gauge-invariant equations. One can see, that a gauge 
function, corresponding to such gauge transformation, should depend on 6 
and z only. After that we introduce the second generating functional 



^p(t) = exp lie f dx p 



c 



dz„a(z, x)F Vfi (z, t)+(3 J dt'(j p (x, t') + r] p (x, t')) 



(20) 

where 3 is a c-number, a(z,x) = \ z x °\ v , C is some fixed closed contour. 
Differentiating $/?(£) by t, applying to it cumulant expansion, differentiating 
once by (3 and putting then (3 equal to (—1), we obtain, due to (5), the two 
last equations of bilocal approximation: 

1 [ y r u e 8 

-<J dz x a(z,y) J dx p a(x,u) — < F Xl y(z,t)F pp (x,t) > - 

\Xq Xq 



y t q 

J dz x a(z,y) J dt'— 



< F x „{z,t)j p {u,t') > + < F Xv (z,t)r)^(u,t') > 



dt 



— J dx p a(x, u) J dt 
xo o 

y 

- J dz x a(z, y) 

X 



< F ppj (x,t)j iy (y,t') > + < F pli (x,t)rj v (y,lf) > 



< F Xu (z,t)j p (u,t) > + < F Xu (z,t)r]^u,t) > 



dx p a(x, u) 



x 



< Fpii(x,t)j v (y,t) > + < F Pfi (x,t)rj v (y,t) > 



+ 



+ J dt' < jv(y,t)j p (u,t') > + < jv{y,t)j»{u,t) > + < j„(y,t)ri P (u,t') > + 



+ < ju(y, t')%(u, t) > + < j^u, t)rj u (y, t') > + < j^u, t')r) y {y, t) > 



+ 



+<W% - u ) = -^-\J dz x a{z, y) < F pli (u, t)F Xv (z, t) > 



du 



P Uo 



/ dt' < F m (u,t)j v (y,t') > + < F m {u,t)r) v {y,t') > , (21) 

A L -J ' 



d 



J dz x a{z,y)— < F Xv {z,t)j p {u,t') > + < F Xv (z,t)r)^u,t') > 

Xq 

- < 3v(y,t)in(u,t') > - < jviy^ri^t') > - 

- < j^{u,t')r] v {y,t) > +26 lu ,6(y - u)S(t - t')- 



e I d 

-— j dv^ <f dw a I J dz x a(z, v) J dx p a(x, w)— < F x ^(z, t)F pa (x, t) > - 



C C 



x 



x 



< F K (z,t)j a (w,V) > + < F x( (z,t)r]„(w,t') > 



W t 

— J dx p a(x, w) J di 
xo o 



dt 



< F pa (x,t)j^(v,t') > + < F p(T (x,t)riz(v,t r ) > 



v 

-J dz x a(z,v) < F x ^(z,t)j a (w,t) > + < F x ^(z,t)ri a (w,t) > 

x L 



dx p a(x, w) 



x 



< F pa (x,t)j^(v,t) > + < F pa (x,t)ri^v,t) > 



+ 



+ / di' < jf(vj)j a (w,t) > + < j^v,t)j a (wj) > + < jz{v,¥)ri a (w,t) > + 



+ <i^>%K?) > + < j a (w,?)riz(v,t) > + 



+ 2(^0 - w) 



- y 

J dz c a(z,y) < F Cl/ (z,t)j p (u,t') > + < F c „(z, fj^u, t') > 

-x a L 



- jdt" < j„(y,t')j^u,t") > + < MvJHW) > + 



+ < MuXftAvJ) > +2M(y - u)5{t" - t') 



_d_ 



< F Xp {u,t)j v {y,t') > + < F Xfl (u,t)r] u (y,t') > 



(22) 



where the terms with space-time derivatives are put to the right hand sides. 
In derivation of equations (21) and (22) we used the following formula [9]: 



<e A B>=<e A > + £ — <C,4 n 5>, (23) 

V n =i n\ J 

where A and B are two arbitrary operators. 

The system of equations (11), (12), (15)-(19), (21), (22) is a complete 

system of equations of bilocal approximation for gauge-invariant correlators 

< ip*(x, f)$(x, x, r)ip(x, t)>, < r/*(x, t?)$(x, x, r)(p(x, t) >, 

< <p*(x, t)$(x, x, t)t}(x, t') >, < F^(x, t)F\p(x, t) >, < F^(x, t)j\(x, t) >, 

< F^(x,t)r]x{x,t) >,< j^(x,t)j v {x,i) >,< j^{x,t)7] v {x,l) > . To obtain 
equations of higher approximations, one should exploit generating functionals 
(10) and (20) in the same way, as it was done for bilocal approximation. 

Hence, there are two alternative ways of investigation of higher correlators. 
E.g., for threelocal correlators, the first way is to derive, using the genera- 
ting functionals (10) and (20), an exact system of equations, where threelocal 
and bilocal cumulants are considered on the same footing. The other way 
(based on the additional assumption, that bilocal cumulants are dominant) 
is to obtain bilocal cumulants from the equations (11), (12), (15)-(19), (21), 
(22) and then to substitute them into the system of equations for threelocal 
cumulants. 

Note, that ANO strings enter equations of bilocal approximation only 
through an additional term in equation (17). Without writing explicitly 
twenty equations of threelocal approximation for correlators < (p*$ip >, 

< rj*<f>ip >, < ^$77 >, < F^Fxp >, < F^rjx >, < F^jx >, < j^j v >, 

< 3 flu >, < tp^tprin >> < tf^Wn >, < <P*^Wlii >, < <P*$<P3n >, 

< rfQipjp >,< ^^ijjp >,< F^Fxpja >,< F^FxpTja >,< Fp V F X pF ai > 

,< jfjujx >,< JuJuVx >,< jiiVvVx >, one can see , tnat m tne correlators 

< Fp V Fx p 3o >, < Fp V Fx p T] a > and < F^F X pF^ > there are groups of terms, 
consisting of bilocal correlators, which have the coefficients, depending on the 
ANO strings. Therefore, ANO strings play nontrivial role in the process of 
formation of the hierarchy of correlators, connecting correlators of different 
orders. 



3 QCD with Spinless Quarks at T = and T > 

In this section we demonstrate, how one can apply the methods of section 
2 to derivation of equations for correlators in QCD with spinless quarks, i.e. 
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the quarks, for which it is possible to omit the term f <7 MZ ,-F MZ , in the quadrized 
Dirac equation, where = §[7^,7^]- Then we generalize this approach to 
the case of finite temperatures. 

The main problem in stochastic quantization of fermions is that the solu- 
tion of naive Langevin equation, obtained directly from the fermionic ac- 
tion, in the physical limit, t — > +00, has an asymptotics e - mt + l P x ^ i. e . 
badly defined in the chiral limit. To avoid this difficulty, the stochastic 
process is modified by introducing a kernel into the Langevin equation: 
1 / dyK(x, y) jfr-p; , where K{x,y) should be chosen in such a way 



as to precisely cancel the negative eigenvalues of j=r A . For example, in [14] 



Stp{x,t) J y v ' ^ J 5^{y,t)' 

Acraf .iirp p>i crprnT-al hps r»f 

6ip(x,t) ' 

K(x,y) = (i-^ + rnj 5{x — y), and the asymptotics of the solution becomes 
e -(p 2 +m 2 )t+t P x ag - n ^ e bosonic case. Our goal is to choose such a kernel, 
that ensures explicit gauge invariance of the right hand sides of equations 
for gauge-invariant vacuum correlators, containing quarks, antiquarks and 
Grassmann stochastic noises. 

In order to get desirable asymptotics of the solutions at t — > +00 (not 
(—00)), i.e. to deal, as usual, with the retarded Green function of the 
Langevin equation, we shall use the Euclidean 7-matrices, multiplied by "i": 

7 = (3a, 74 = i(3, where a = ^ ^ j , (3 = ^ J ^ j , (24) 

which satisfy anticommutational relations {7^,7^} = —25^ v . 

Let us introduce the kernel K(x,y,t) = (iD(x,t) + m)6(x — y) into the 
naive Langevin equation, obtained from the Euclidean QCD action with 7- 
matrices, defined in (24), 



S = J dx 



where D fl = d, - igA, 



For spinless quarks one gets the following Langevin equations: 

ip = {D 2 -m 2 )ip + e, (25) 
ij, = (V 2 - m 2 )i) + 0, (26) 

where ip = ^ + 74, T>^ = + igA^. The correlator of Grassmann noises 
changes correspondingly: 

< 9 a (x,t)9^y,t') > 6 o =2[ib(x,t)6(x-y)+m6(x-y)] aP 6(t-t l ), (27) 
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where for an arbitrary functional of the 9's we have the stochastic expectation 
values: 



<F[e,e\> 6 - e = 



J D9D9F[9, 9]exp[-\j dxdt9(x, t)(W + m) _1 0(a;, t)] 



J D9D9exp[-± I dxdt9(x, t) (iD + m)- l 9(x, t)] 

(28) 

Solving (25) and (26) via Feynman-Schwinger path integral representation 
with the initial conditions ip(x, 0) = ip(x, 0) = and using (27) and (23), one 
obtains in bilocal approximation: 

< 1pp(x,t)$(x,X,T)lp a (x,t) > = 



min(t,t) 

= 2 / dt' f dy(Dz),yK- z (i,t') 
o 



^{Dz) xy K z {t,t') 



+ 



y=y 



+ m(Dz) xy K z (t,t% = v } < Wy xx y >, 



(29) 



where < ... > means the averaging over the vacua of rf^ and (9,9),!]® is a 
Gaussian stochastic noise with the bilocal correlator < rj^(x,t)rjl(y,^) > v = 
= 25^5 ab 5{x - y)6{t - t'). Here 



$>(x, x,t) = P exp 



X 

ig J dw^A^(w,T) 





r y 


< ig 


/ 




L x 



d^(z(0,0 + / dw^(w,T) + / dz^(z(0,0 



Wyxxy = Pexp { 

(30) 

To perform differentiation in the equation (29), let us extract explicitly 
the dependence on the point y from the measure of integration and from 
the kinematical factor K z (t,t'). To this end we pass to the integration over 
trajectories 

y-x 



u{0 = m + j-^tt-t')-y. 



(31) 



After that differentiation may be easily performed, and the result has the 
form 
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min(t,t) 

2 / di' jdy(Dz)- X yK- z (t,t') 
o 



i x — y ( x -v) 2 
- e i(f-t) 

2t' -t 



Du) 



00" 



K u {t,t') < Wyx X y > +mj {Dz) xy K z (t,t , ) y= y < Wj 



yxxy 



> 



a[3 



(32) 



Here in the first Wilson loop in the right hand side it is implied, that in the 
last term in (30) z M is expressed through via (31). 

In analogous way one may obtain equations for two other gauge-invariant 
correlators, containing spinor fields, in bilocal approximation: 



< ipp(x, t)$(x, x, r)6 a (x, t) > and < 0@(x, t)$(x, x, r)ip a (x, t) > . (33) 

In contrast to the AHM, in the case of QCD it is, of course, impossible to 
calculate exact Green function of the Langevin equation for F£ v (x,t) (this 
problem is equivalent to the solution of Yang-Mills equation in the general 
form). Therefore, in order to obtain all the rest equations for gauge- invariant 
correlators in bilocal approximation, one should act in the same way as it was 
done in section 2 for derivation of equations (21) and (22). Using Schwinger 
gauge and introducing the generating functional 



<3>/3 = Pexp < 



ig jdx^ j dz v a{z,x)F V[l {z,t) + (3 f dt '(gj»(x,l/) +rj fi (x,t')) 

k C \ 

(34) 

where j° = ipj^ip is a quark current, differentiating (34) by t and applying 
to it cumulant expansion, one obtains, due to the Langevin equation 

A; = (D x F x ,y + 9J ; + n% 

where D\ is the ajoint covariant derivative, the first equation. After that one 
should differentiate it five times by (3, use the formula (23) and in all the 
equations put (3 equal to (—1). Finally, to obtain a minimal closed system of 
equations of bilocal approximation, one should avoid D\ for d\ in the right 
hand sides of these equations. This may be done, using the expressions for 
higher orders' path-ordered cumulants [3,10] and the formula 
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tr <<f>(x ,u,t)(D x F Xfl (u,t'))<!>( 

\ d } 

= tr < T7— < F Xfl (u, x , t)G Ml ... Mn > +ig / dx p a(x, u)- 

< F Xfl (u, x ,t)G ^..^ n F Xp (x, x ,t) > - 

- < Fx ll (u,XQ,t)Fxp(x,XQ,t)G tll .. mlln > j, 

where F^x, Xo, t) = <$>(xo,x,t)F fJ , u (x,t)<$>(x,XQ,t), any parallel transporter 
between xo and any other point goes along straight line, G^...^ = 
= C Ml ... Mn (xi, ^i, x n , t n , xq, £) is, generally speaking, a product of some num- 
ber of F^j^ and rj^, which are given in the points at the moments 
ti, ...,t n of Langevin time respectively, and all the parallel transporters 
between xo and each of these points are given at the same moment t. This 
method of derivation of the system of equations in the case of pure gluody- 
namics was used in [10]. 

After that, substituting < F fJ , l ,(x,Xo,t)F\ p (y,xo,t) >, obtained from this 
system of six equations for correlators 

< F^(x, x , t)F Xp (y, x , t)>, < F^x, x , t)j\{y, x , t, t') >, 



< F^(x, x , t)r}\(y, x , t, t') >, < j M (x, x , t, t')j„(y, x , t, t") >, 

< j»(x, x , t, t')r) u (y, x , t, t") >, < ^(x, x , t, t')r] v {y, x , t, t") >, (35) 
where 

j M (x, x , t, t') = $(x , x, f)j M (x, ^o, 

r^(x, x , t, t') = $(x , x, f)r/ M (x, 0^(^ 5 ^o, t), 

into the averaged Wilson loops, expanded up to the second cumulant via 
nonabelian Stokes theorem in (32) and in the right hand sides of equations 
for correlators (33), we obtain all the correlators of bilocal approximation. 
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In order to obtain equations for higher correlators, one should introduce 
the analog of generating functional (10) and follow the algorithm, described 
in section 2. 

To conclude this section, we generalize the method of derivation of the sys- 
tem of exact equations, suggested above, to the case of QCD at finite temper- 
atures. It was shown in [15], that the ordinary results of finite temperature 
field theory may be reproduced in the framework of stochastic quantization, 
if one requires the stochastic noise to be periodic for bosons (antiperiodic for 
fermions) in the coordinate £4 with the period f3 = ^. Therefore, we supply 
equations (25) and (26) with the following conditions: 

6(x, X4, t) = — 6(x, X4 + (3, t), 6(x, X4, t) = — 6(x, £4 + (3, t) (36) 

and modify the meaning of the averaging over (9, 9) from (28) to 



<F[9A >ee= 



J D9D9F[9,9]exp 


r 1 p 1 i 

-Udx 4 J dxdt9(x, t) (iD + m)- 1 9(x, t) 



J D9D9exp 


- p _ 

— 2 / dx 4 J dxdt9(x, t)(iD + m)~ 1 9(x, t) 





Then it is easy to check, that in the absence of gluonic fields, the solution of 
equation (25) with initial condition ifj(x, 0) = has the form 

00 P 

if)(x,t) = J dt'J dy A J dyG(x-y,t-t')6(y,t'), (37) 



where 



0(t) +2? , dk j, x _, kl 



ik n x—(k 2 l +m 2 )t 



(3 n^ooJ (27T) 3 
* ft ~\ - ( 2n + 1 ) 7r 

k n = {k,u n ), u n = . 

Solving equation (26) for a free antiquark with the initial condition ip(x, 0) = 
= in analogous way, one gets after some calculations 



- 1 r dk * 1 

< ip p (x',t')ip a (x,t) > e - e = - E / TTT*{-kn + m) a pe lkn{x ~ x,) - 

P n=-oo J \^) 6 
e -{kl+m 2 )\t-t'\ _ e -(kl+m 2 )(t+t') 

kl + m 2 



15 



that in the physical limit t = t' — > +00 yields standard result for the propa- 
gator of a free fermion at finite temperature [16]. 

In the case, when gluonic fields are present, one, solving equation (25) with 
the initial condition ip(x,Q) = through Feynman-Schwinger path integral 
representation, obtains instead of (37): 

t 

il>(x,t) =jdt'j dy(Dz)™ y K z (t,t')<S>(x,y)6(y,t), 


where 



( Dz )xy = lim n 77=^ E -rAf Wl 

y n^oo l=1 (4:7re) z n= _oo (27r) 4 



C„(fc) = z,(ek) - z,(e{h - 1)) [17], 



<$>(x, y) — Pexp 



ig / dzpA^zte),?) 



Therefore, one can check, that, for example, equation (32) changes as follows: 



i(t,t) 



< 



ipp(x, i)$(x, x, r)t/j a (x, t)>=2 J dt' J dy(Dz)^ y K z (t, t') 



i x — y (jE-vr 

- e 4(F^t) 

2f-t 



■ J (Du) w K u (t, if) < W rxX y > +m / {Dz) w xy K z {t, < W yxxy > , 

u(tf)=u(t)=0 af3 

(38) 

where 

(Du) w = lim n jrH E W / , (39) 

A(fc) = u(efc) - u(e(fc- 1)). 

Solving equations for correlators (35), one should keep in mind, that the 
stochastic noise fields r}^(x, t) are periodic in the coordinate £4 with the period 
P. 
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4 QED: Including Spin Effects 



In this section we shall demonstrate, how one can exactly take into account 
spin effects in the framework of suggested approach. Let us consider for 
simplicity the case of QED and, instead of equation (25), solve exact Langevin 
equation 



where 



^ = _(i) 2 + m 2 )^ + 5 



D 2 = -(D 2 + -o, v F, v ), 



(40) 



via Feynman-Schwinger path integral representation. We shall follow the 
method, suggested in [11,12] for derivation of a propagator of a spinning par- 
ticle in the external Abelian field. The retarded Green function of equation 
(40) has the form: 



G(x, y, t) = 9{t) j (Dz) xy Dp exp lijdr 



pz + i(pu - eAJz,r)) 2 + 



+-F ltv (z,r) 



+ im A 



| Texp J 



P\{t')i\(It' 



J \p=o 



where p^ are Grassmann sources, anticommuting with 7-matrices. Using 
Wick theorem [18], we get: 



Texp 



■ L 

j P\(t')i\<It' 



= exp 



t t 



- J dn J dr 2 p^{Ti)sgn(n - T 2 )p M (r 2 ) exp J p\{r)j\dT 







/ D£ exp \ f dr 



J Dt; exp 



\ldrUr)Ur) 



■exp^^— Jexp 



-% J p fl (r)9^dT 



0=0 
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where £ M (r) are odd trajectories. To simplify this expression, one can pass 
to the integration over trajectories x(r) = |(£(t) + 0) and obtain: 



G(x,y,t) = 9{t) j (Dz) xy expli-f fi 



_d_ 

W, 



J Dx exp \Jdr 



Xm(«)+Xm(*)=^ 



2: 



m + 



\0u=O 



In order to avoid the restriction to the trajectories, over which we perform the 

integration, let us pass to the integration over trajectories u(r) , such that 
t 

x(r) = I / sgn(r — r')u(r')dr' + |6*. Noticing, that for an arbitrary function, 
given at the Grassmann algebra, the following formula takes place: 

( ' 9 )f{°) \e=o= f(^) ex v{K^l^\^ 



expiry p 



06 



/// 



one obtains the following expression for G(x,y,t): 



( t 



G(x,y,t) = -0(t) J {Dz) xy exp< J dr 

{0 



z 



m 2 + ieA^Zjrfzp 



le 



-—F^(z,t) 





' t 


| j Du exp < 


jdr 




lo 


t 





le 



■sgn{T - t')u^t')F^{z,t) J dr" sgn{r - t")u v (t") + 



+ 9 / dT ' s 9 n i T ~ t')u^{t')u^{t) - y j dTsgn{r - t')u^{t')F^ v {z,t 



d 



dCu 



26 \ 1 



-- j dr'] dr"u,{T')sgn{r' - t")u,(t") 


/ Du... may be represented in the form 



>exp« M 7M)k=o- ( 41 ) 



Du exp 



^u^T^z, e)u v + I^u^ 
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where 



3 c) 1 d 

Tpviz, e) = ieeF^ v (z)e + -S^e, 1^ = ieF^(z)e— + 

and we used the following notation: for example, 

t t 



coeu = J dr J dT"u){T')sgn{r' — t")u{t"). 







Here and everywhere later, where such a notation is used, the Langevin time 
argument of is omitted, since it is integrated over. Then 



detT(z, e) 



1/2 



detT(z,0) 
= exp I i J de'tr[Q^(z, e')F^(z)] 



exp^—I^T \z,e)]^I^ = 



exp 



where 



Q^(z,e) - e[T l (z, e)]^e, R^{z,e) - [T 1 (z,e)] fiV . 
Differentiation in (41) with the help of the formula 

d 2 



exp ( B^y 



-j exp(i(\j\)\ (=0 = 1 - iu^B^ + i^B^B*^, 



where 

75 = —271727374, "*" denotes the dual tensor, 
yields the final result: 



G{x,y,t) = —expl- 



e(e - 2i) 



1 



J (Dz) xy exp < J dr 



— 171 + 



+ieA fi (z, t)z^ - ^a^F^(z, r) + ^75 J dr'F^(z, t)F*„(z, t') 



+ 



+1 J de'tr[Q^(z, e')F^(z)] + A{z, e) 



where 
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A(z,e) = - l -e 2 



a Xp F pX (z)eR pi/ (z, e)F vp (z)e - j 5 (F^ x (z)eR^(z, e)F vp (z)e)- 



■{F aX (z)eR af3 {z, e)Fp p (z)ey 



e 
4 



a Xl ,F^ x (z)eR pi/ (z, e) + o w R pv {z, e)F vp {z)e- 



-7fa(Fii\(z)eRiiv(z, e))(F aX (z)eR av (z, e))* - ^{R p „{z, e)F vp (z)e)- 



{R pa {z,e)F ap (z)e) 



+ 



8 



o pv R pv {z, e) - 75^(2:, e)R*^{z, e) 



In the same way one may obtain the Green function of the Langevin 
equation 

^ = _^ (t3 2 + m 2 ) 

where 

To generalize the equation (32) to the spinor case, one should again, for 
the group of terms in the functional integral, pass to the integration over 
trajectories (31). The result has the form: 

min(t,t) 

< ipp(x,t)$(x, x, T)ip a (x, t) >= 2 J dt' j dy(Dz)xyKz(t, t', -e)- 



^—pr(x - y) l5 e^ J (Du) 00 K u (t, t', e)- 



2{f-t) 



< W rxX y \ex V r-Jd^ l-a^F^z, + 75 / £)F;„(z, £') 



+ 



+ % j de'tr[Q pi/ (z, e')F pv {z)} + A(z, e) 



K?-t')+«(«)+y 



exp \ - j d£ 
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+ 



M J de"tr[Q Xp (z, e")F Xp {z)} + A(z, -e) ^ > + 

5/3 



+m y (Dz) x gK z (t, t\ e)| y= _ < j ex]? j 1 1 d£ 



c -1 o 

+75 / + * / de'trlQ^z, e')F^(z)] + A(z, e) 

v J 



0:7 



exp \ - J d£ 



*XpF\p(z,£) - 75/ deF Xp (z,OF* Xp (z,e) 
t> 



+ 



+1 J de"tr[Q Xp (z, e")F Xp (z)] + A(z, -e) 



} > 

7/3 J 



(42) 



where 



K z (t,t',e) = ^exp 



'z 2 



J -dt-m\t-t') + --ie- 
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Therefore, we reduced the right hand side of the equation for 
< iljp(x,t)&(x, x,T)ifj a (x,t) > to the path integral over bosonic trajectories 
from the functional, whose spinor structure is given explicitly. The following 
steps of derivation of the system of equations for gauge-invariant correlators 
are similar to the previous two sections and based on the introduction of the 
corresponding generating functionals, analogous to (10) and (20). 



5 Conclusion 



In this paper, using stochastic quantization method [1], we developed a 
systematical procedure of derivation of the systems of equations for gauge- 
invariant vacuum correlators in different field theories. In section 2 on the 
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example of AHM we demonstrated, how this method may be applied to 
quantization of classical solutions (ANO strings); the quantization of more 
complicated objects will be treated elsewhere. The generating functionals 
(10) and (20), introduced in this section, appropriately modified, may be 
used to derivation of equations in any other gauge theory with the help of 
the algorithm, described in sections 2 and 3. This algorithm, was used in 
section 2, where the minimal set of equations of bilocal approximation was 
derived, and the role of ANO strings in the higher correlators was discussed. 
The suggested approach yields the gauge-invariant method of quantization 
of the theories with spontaneous symmetry breaking; this problem will be 
investigated in the future publications. 

According to the general principa of the Method of Vacuum Correlators, 
the obtained system of equations is splitted into two parts: the part, con- 
taining pure gauge fields' correlators and matter fields' currents and the part, 
containing matter Green functions, so that the equations of the first subsys- 
tem may be solved independently from the second equations, but not vice 
versa. The connection of this fact with the large- A" behaviour of QCD and a 
new type of the Master field equation, based on the suggested approach, will 
be treated elsewhere. 

The method of derivation of equations in QCD with spinless quarks was 
generalized at the end of section 3 to the case of finite temperatures, and in 
section 4 we showed on the example of QED, how one can take into account 
spin effects exactly. 

The equations, obtained in the case of QCD, contain both perturbative 
and nonperturbative gluonic contributions to the vacuum correlators. The 
problem of separation of perturbative gluonic contributions in all the terms 
of cumulant expansion and the regularization of the obtained equations will 
be the topics of separate publications. 

The results, presented in this paper, were partially reported at the In- 
ternational Workshop "Nonperturbative Approaches to QCD", Trento, Italy, 
July 10-29, 1995. 
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Appendix 

Derivation of the Green function (8) 

Below we present explicit derivation of a Green function of the Langevin 
equation 

Fnv{x,t) +d v d\Fx li (x,t) - d^d x F Xv {x,t) = f^ v (x,t), (A.l) 

where 

f^(x,t) = d^j^x.t) + rj u (x,t)) - d v (jp(x,t) +rj^x,t)). 

Writing equation (A.l) in the momentum representation: 

F^(k, t) + k^k x F Xv (k, t) + k v k x F^ x (k, t) = f^(k, t), (A.2) 

one can check, that the retarded Green function G Xpia /3(k, t) of equation (A.2) 
satisfies the following equation: 



Q 2 

^-^,\ P -q^ + k (1 — P) M jy,Ap 

Here we introduced the projectors 



G\p,ap(k,t) = l^, a pS(t). (A.3) 



ljuz^Ap = ~^{^ii\$vp ~ 5 pp 5 vX j and P^,Ap — — ^T^xT^p — Tp P T uX j, 

where 

rp c _ kpk v 

J-pv — Ofu, 2 , 

possessing the following properties: 

~^-pv,\p ~^-i/fi,\p ~^-pv,p\ ^-Xp^pvt ~^-pv,\p~^-\p,a>l3 ~^-p,v,aj3 (-^-•^O 

(the same properties hold true for P^Ap), 

P^,Ap(l - P)Ap,a/3 = 0. 

Using (A. 4) and (A. 5), one can easily get G^ a /3(k,t) from equation (A.3): 
G^ a(3 (k,t) = 0(t)[(l - P)e- pt + P] w/3 . 
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Returning back to the coordinate representation, one obtains, that the 
lution of equation (A.l) with the initial condition A^(x,0) = A^(x) has 
form (7), where 



= ^-J dk e^k c 
that after some calculations yields (8). 



(l-P)e- fc * + P 
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